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The Miiller-Israel-Stewart second order theory of relativistic imperfect fluids based on Grad's mo- 
ment method is used to study the expansion of hot matter produced in ultra-relativistic heavy ion 
collisions. The temperature evolution is investigated in the framework of the Bjorken boost-invariant 
scaling limit. The results of these second-order theories are compared to those of first-order theories 
due to Eckart and to Landau and Lifshitz and those of zeroth order (perfect fiuid) due to Euler. 
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High energy heavy-ion coUisions offer the opportunity 
to study the properties of hot and dense matter. To 
do so we must follow its space-time evolution, which is 
affected not only by the equation of state but also by 
dissipative, non-equilibrium processes. Thus we need to 
know the transport coefficients such as viscosity, thermal 
conductivity, and diffusion. We also need to know the 
relaxation coefficients. Knowledge of the various time 
and length scales is of central importance to help decide 
whether to apply fluid dynamics or cascade or a com- 
bination of the two. The use of fluid dynamics as one 
of the approaches in modeling the dynamic evolution of 
nuclear collisions has been successful in describing many 
of the observables [|l],|j . So far most work has focused on 
the ideal or perfect fluid and/or multi-fluid dynamics. In 
this work we apply the relativistic dissipative fluid dy- 
namical approach. It is known even from non-relativistic 
studies 1^ that dissipation might affect the observables. 
The first theories of relativistic dissipative fluid dynam- 
ics are due to Eckart Q and to Landau and Lifshitz 
These theories are now known to have some undesirable 
features: they lead to Navier-Stokes-Fourier laws which 
are parabolic in structure and therefore may propagate 
signals with speeds exceeding that of light. A qualitative 
study of relativistic dissipative fluids for applications to 
relativistic heavy ions collisions has been done using these 
first-order theories 

Second-order theories of dissipative fluids due to Grad 
PI , Miiller [Q, and Israel and Stewart were intro- 
duced to remedy some of these undesirable features. In 
second-order theories the space of thermodynamic quan- 
tities is expanded to include the dissipative quantities for 
the particular system under consideration. These dissi- 
pative quantities are treated as thermodynamic variables 
in their own right. The phenomenological formulation of 
the transport equations for the first-order and second- 
order theories is accomplished by combining the conser- 
vation of energy-momentum and particle number with 
the Gibbs equation. One then obtains an expression for 
the entropy 4-current, and its divergence leads to entropy 
production. Because of the enlargement of the space of 
variables the expressions for the energy-momentum ten- 
sor T^'^, particle 4-current TV^, entropy 4-current S^, 



and the Gibbs equation contain extra terms. Transport 
equations for dissipative fluxes are obtained by impos- 
ing the second law of thermodynamics, that is, the prin- 
ciple of nondecreasing entropy. The difference between 
the two stems from the entropy 4-current: the standard 
irreversible thermodynamics of Eckart-Landau assumes 
that the entropy 4-current should include terms linear in 
dissipative fluxes and hence they are referred to as first- 
order theories of dissipative fluids. On the other hand the 
extended irreversible thermodynamics of Miiller-Israel- 
Stewart includes terms quadratic in dissipative fluxes and 
hence they are referred to as second-order theories of dis- 
sipative fluids. The kinetic approach is based on Grad's 
14-moment method Jl]| . The resulting equations are hy- 
perbolic and lead to causal propagations ||l^,|l^. For a 
review on generalization of the Miiller-Israel-Stewart the- 
ory to a mixture of several particle species see p^ . 

The formulation of relativistic hydrodynamics can be 
found in standard text books • The energy equa- 

tion governing fluid motion is given by 

De^-{e+p + n)e TT^'" v^' w/x - V/.9^ + Q^a^ , (1) 
where 



(energy density) , 
(pressure) , 

(convective time derivative) , 
V*" = A'^'^di, ( gradient operator) 

daU°' (volume expansion ) , (2) 
(4-velocity) , 
Du^ (4-acceleration) , 
(viscous pressure) , 
(viscous stress tensor) , 
q'^ (heat flow) , 

A^'' = gf^"^ — u^u'^ is the projection tensor orthogonal 
to the 4-velocity and g^'^ = diag(+l, —1, —1, —1) is the 
metric tensor in Minkowski space-time. In the Eckart 
theory the dissipative contribution to the bulk pressure, 
the heat flux and the shear viscous tensor are given by 
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where C , , A are the bulk, shear and thermal conduc- 
tivity coefficients respectively. They are required to be 
positive by the second law of thermodynamics. T is the 
temperature of the system. The angular bracket notation 
is defined by: 



i (A;^ + A^! ^l) - i A^"- A,, 



A"^ . (6) 



The above expressions for dissipative fluxes can also 
be obtained in kinetic theory by employing the first 
Chapman-Enskog approximation fl^ . In second-order 
theories we have to solve for these dissipative fluxes from 
their evolution equations. One can show [p^ that in 
the 1-1-1 Bjorken scaling solution hypothesis ||l6|, the 
system of transport equations given by fl^ j becomes 
tractable. The same simplifications occur in 1-1-1 dimen- 
sional Bjorken hypothesis in first order dissipative fluid 
calculations In 2-1-1 and 34-1 dimensional flow these 
equations become much more complicated. This is a sub- 
ject of current study fl^ . 

I give here the second-order equations for 1-1-1 dimen- 
sional scaling solution. These equations were derived 
from kinetic theory using Grad's 14-moment approxima- 
tion method ||ll[]. We will only need the Miiller-Israel- 
Stewart [13| equations in the following form: 



where 



To'' 



1 /v^r 



f3i\ T 



(7) 

(8) 
(9) 



1 

T2 



(10) 



are the relaxation times, /3o , Pi and /32 are the relaxation 
coefficients. These three new coefficients are functions 
of primary thermodynamic variables like pressure, num- 
ber density and energy density and hence depend on the 
equation of state. Relaxation time is the distinguishing 
feature of hyperbolic causal theories which is not present 
in the first-order theories. Here r, is the time taken by 



the corresponding dissipative flux to relax to its steady- 
state value. 

For the 1+1 dimensional Bjorken | [l6t similarity fluid 
flow the energy equation becomes: 



de 



P 2 1 
$ 

3t 



n- = , 

T 



(11) 



where $ = tt — tt^^ is determined from the shear viscous 
tensor evolution equation 



d 1 1/1 1 d /f32 
— $ = $ _ _$ _ ^ T — — 

dr 2 \t f32 dr \T 



The equation for 11 will not needed as explained below. 
For this initial study a simple equation of state is used, 
namely that of a weakly interacting plasma of mass- 
less u ,d ,s quarks and gluons. The pressure is given 
by p — e/3 = aT^ with zero baryon chemical poten- 
tial. Here a is a constant determined by the number of 
quark flavors and the number of gluon colors. In the 
case of massless particles the bulk pressure equation (0) 
does not contribute since the bulk viscosity is negligible 
or vanishes Q . The only relaxation coefficient we need is 
P2 which for massless particles is given by P2 — 3/(4p). 
The shear viscosity is given by 77 — bT^ where 6 is a 
constant determined by the number of quark flavors and 
the number of gluon colors. The energy equation (|^) and 
the viscous stress tensor equation (|o|) can be written as 
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Here Nf is the number of quark flavors, taken to be 3, 
and as is the strong fine structure constant taken to be 
0.5. For a perfect fluid and a ffi'st-order theory the energy 
equation (Ww can be solved analytically to give 



T(r) = To 
T{t) = To 



1/3 
1/3 



1 



(perfect fluid) , 



6 a To To 
(first-order theory) , 



2/3 



(17) 
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In the first order theory we do not have the relaxation 
coefficients. Then equation (|l2|) gives $ (4?7/3)/r. 
The above equations and the numerical solution to the 
second-order equations ( p^ ) and ( p^ are presented in 
Figs. |l| through m We choose the initial temperatures to 



2 



correspond to those expected at RHIC and LHC, namely 
To=500 MeV at RHIC and ro=1000 MeV at LHC. In 
Figs. |l| and || the initial time tq is estimated by using the 
uncertainty principle |^^: tq ■ {E)q ~ 1 where {E)^ ^ 3To 
for massless particles. This results in to=0.13 fm/c at 
RHIC and to=0.07 fm/c at LHC. The initial value used 
for <i>, which must be specified independently for the sec- 
ond order theory, is taken to be p/3. We choose this value 
since the second order theory is based on the assumption 
that the dissipative fluxes are small compared to the pri- 
mary thermodynamic variables, namely p, n,and e. How- 
ever, a thorough study of the initial conditions on these 
dissipative fluxes is needed and should perhaps be found 
from microscopic models. This is as subject of current 
study jl^ . The effect of dissipation is more pronounced 
at the very early stages of heavy ion collisions when gra- 
dients of temperature, velocity, etc., are large. At late 
times the effect of dissipation vanishes. In Figs. ^ and 
^ we take a constant initial time ro=1.0 fm/c which is 
the characteristic hadronic time scale. Euler hydrody- 
namics predicts the fastest cooling. The first-order the- 
ory significantly under-predicts the work done during the 
expansion relative to the Miiller-Israel-Stewart and Euler 
predictions. Thus the temperature decreases more slowly 
with the inclusion of dissipative effects. This would lead 
to greater yields of photons and dileptons, and also the 
transverse energy /momentum would be reduced as the 
collective velocities are dissipated into heat. The system 
takes longer to cool down. This will delay freeze-out. 
Also entropy, s — AaT'^, is enhanced. This is important 
because entropy production can be related to final mul- 
tiplicity. With respect to entropy production due to par- 
ticle production, see [2l] ]. Given some initial conditions 
we want to investigate the importance of second order 
theories as compared to first order theories and perfect 
fluids. Let us now analyze the differences between the 
second-order and first-order theories. The first thing we 
notice is that the Eckart-Landau theory predicts that at 
early times the temperature will first rise before falling 
off. This is more pronounced when we have small ini- 
tial times. Naively one would expect that the system 
would cool monotonically as it expands even in the case 
of dissipation where energy-momentum is conserved. On 
the other hand, it is seen that for large initial times and 
high temperatures the three theories have a similar time 
evolution. As can be seen from Fig. ^, all three cases 
start at the same point and then fall off with time. The 
difference stems from the fact that in the second-order 
theory the transport equations of the dissipative fluxes 
describe the evolution of these fluxes from an arbitrary 
initial state to an equilibrium state. The first-order the- 
ory, though, is just related to the thermodynamic forces 
which, if switched off, do not demonstrate relaxation. 
Hence they are sometimes referred to as quasi- stationary 
theories. As can be seen from Fig. ^ it is before the 
establishment of an equilibrium state that the two the- 



ories differ significantly. In ultra-relativistic heavy ion 
collisions, where the fluid evolution occurs in very short 
times, the second-order theories should be used to ana- 
lyze collision dynamics. In doing so a full analysis re- 
quires that all the dynamical equations with more real- 
istic equations of state and transport coefficients be con- 
sidered 

I will conclude by pointing out some of the advantages 
and challenges of the second-order theories. Second-order 
theories, being hyperbolic in structure, lead to well-posed 
initial- value (Cauchy) problems. They also lead to causal 
propagation. Unlike the first-order theories, second-order 
theories have relaxation terms which permit us to study 
the evolution of the dissipative fluxes. The challenge we 
face is the increase in the space of thermodynamic vari- 
ables. We now have, in addition to the transport coef- 
ficients, new coefficients in the problem. These are the 
relaxation coefficients /3i and the coupling coefficients a^. 
These new coefficients depend on the primary thermo- 
dynamic variables, such as n, e and p and therefore are 
determined by the equation of state. Like viscosity and 
thermal conductivity, which are required to be positive 
by the second law of thermodynamics, these new coeffi- 
cients are constrained by hyperbolicity requirements. In 
principle , in order to solve the second-order relativis- 
tic dissipative fluid dynamic problem, one still needs the 
equation of state, initial conditions and the transport co- 
efficients. 

To probe non-equilibrium properties of matter pro- 
duced in heavy ion collisions we need a non-equilibrium 
fluid dynamics model to analyze observables. The rela- 
tivistic fluid dynamics modeling of heavy ion collisions 
will have to include dissipation and thermal conduction. 
One will then have to use the hyperbolic theories of rel- 
ativistic dissipative fluids because of their universality. 
Hyperbolic theories might prove to be convenient in con- 
structing hydro-molecular dynamic schemes [p2| in which 
a phenomenological fluid dynamics model is coupled to 
a microscopic kinetic model such that microscopic kine- 
matic quantities may be obtained. Dissipation mecha- 
nism might be important if we deal with the event-by- 
event based hydrodynamics p^ . 
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FIG. 1. Proper time evolution of temperature for a RHIC 
scenario: To = 0.13 fm/c and To = 500 MeV . 
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FIG. 2. Proper time evolution of temperature for a LHC 
scenario: To = 0.07 fm/c and To = 1000 MeV . 
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FIG. 3. Proper time evolution of temperature for a RHIC 
scenario: To = 1.0 fm/c and To = 500 MeV. 
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FIG. 4. Proper time evolution of temperature for a LHC 
scenario: To = 1.0 fm/c and To = 1000 MeV. 



5 



